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. -. ■ The time-reversal mirror (TRM) of acoustic waves is used as a focusing technique that is robust 

against disorder and chaos. The propagation of an ultrasonic pulse from the focal point is detected 

and recorded by a surrounding array of transducers that enclose the cavity region. The registered 
C3 ! 
'-^ ■ signals are then played-back in the reverse temporal sequence and an imperfect Loschmidt echo 

is obtained. A similar technique developed in the quantum domain, the perfect inverse filter 
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(PIF), ensures an exact reversion. Applied to the acoustical case, it allows an improvement in the 
design of a minimally invasive inverse filter. Although that version of PIF accounts for external 



source condition, now we propose a generalization where the initial excitation is created inside the 
cavity. We present the theoretical principles for this form of PIF and perform a numerical test in 



a multichannel system that improves the focalization quality compared to the TRM version. 



I. INTRODUCTION 

In acoustic experiments, the time-reversal mirror (TRM) technique has been developed 
by Mathias Fink and collab. [1] in order to focus ultrasonic waves through heterogeneous 
media. The remarkable robustness of this technique against inhomogeneities allows novel 
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applications in medical physics [2j, oceanography [3| and telecommunications [4]. The acous- 
tic signals diverging from a localized source point or scatterer are recorded by an array of 
transducers until all excitations have decayed. Their registries are then played-back follow- 
ing the inverted time sequence. Notably, with this simple prescription, the initial excitation 
is recovered in the source region up to a good degree of fidelity. This constitutes an al- 
ternative form of achieving the Loschmidt Echo 5|, i.e. the time reversal of an excitation 
that has spread in a complex medium. In its closed version, called time-reversal cavity, 
the array of transducers B encloses the source region separating the space in two portions: 
the closed cavity A, where the focusing takes place; and the outer region, towards which 
all the excitations finally escapes. In order to ensure the perfect refocusing and quantify 
the stability of the TRM, we have introduced a formal prescription for quantum systems 
denoted as the perfect inverse filter (PIF) [6J. This achieves exact focusing by injecting 
a wave function that has to be corrected to account for the time reversed signal and the 
system response. The injection prescription was derived as a solution to the Schrodinger 
equation with boundary source conditions that arises when a wave packet arrives at the 
cavity from the outside. As the wave packet enters at the cavity, it interferes with its re- 
flected component and finally spreads out in space and time. The reversed counterpart is 
obtained via the injection process and a burst of wave probability is emitted from A. This 
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is the quantum version of the sound bazooka, suggested by M. Fink and collaborators [7]. 
furthermore, we implemented this prescription in a classical chain of harmonic oscillators 
8|. In that case, a low frequency mode (e.g. a heavy mass) was coupled to a semi-infinite 
chain of lighter masses leading to broad continuum spectrum of excitations that acts as an 
"environment". The energy stored in the surface oscillator decays almost exponentially and 
its propagation was registered by a single transducer inside the environment. In this case 
where all the normal modes participate in the local excitation, the PIF correction becomes 
effective and both the dynamical build-up and the final contrast, spatial distribution of the 
recovered initial condition, are substantially improved. 



In addition to the TRM, Fink and collab. proposed the spatio-temporal inverse filter 
(STIF) for acoustic waves as a focusing method that optimizes spatial contrast |9|. It 
requires the inversion of the Fourier transform of the direct wavefield propagator. This 
matrix is build from the temporal response at the transducers to a Dirac delta function in 
each control point within the focal region. Unlike the TRM, the STIF injection acts as a 
highly invasive technique because it is necessary to measure the full propagator whereas the 
former only need to compute the row related to the focal point. Recently, an alternative 
proposal by Vignon et al. 10|, denoted minimally invasive STIF (miSTIF), achieves the 
same focusing quality than STIF. The idea behind this procedure is that all the information 
in the wavefield propagator can be deduced from the backscattering signals of the transducer 
array at the boundaries. Hence the technique is no more invasive than TRM. In this context, 
miSTIF and PIF result similar because they correct the TRM injection function using the 
propagator in the reduced basis of the transducers. 

In this work, we explore the PIF procedure for the case in which the source is placed 
inside the cavity. In Sec. II the general PIF formula is deduced from the properties of the 
system Green's function. The obtained result is then applied in a semi-infinite chain with a 
single transducer and also to a general multichannel system. For this last case, a numerical 
simulation is performed yieding a favourable test of the most difficult task: the time reversal 
of an initially localized excitation. In Sec. Ill the PIF is reobtained from the optical theorem 
and compared with the miSTIF technique. 

II. TIME-REVERSAL WITH INTERNAL SOURCE 

We look for the situation in which at the initial time all excitations are located inside 
the cavity. If we denote by ip x {to) the coefficient of the initial condition at site x, then we 
have 

J(to)=$A(t )=J2Mto)\x)- (1) 

x&A 

During the propagation from t = until a registration time £r, the wave function is detected 
and recorded as crosses the boundaries B 

f B (t) = J]Vs(*) \x) = ihG BA (t)M0), 0<t<t R , (2) 

xeB 



where <&B,A{t) is a matrix arrangement whose element G(x, t; x', 0) relates the response at 
site x £ B to a Dirac delta function <5(£) in i' G A following the Schrodinger equation whose 
Hamiltonian is H 

ih^G(x, t; x', 0) - Y, H x , x „G(x", t; x', 0) = 5 XjX >5(t). (3) 

The Fourier transform of this matrix gives the retarded Green's function 

G(x,x , ;s)=[(eI-m)- 1 ] XtX ,. (4) 

Assuming that there are no localized states, the registration time t R is taken long enough 
in order to ensure that all excitations have been left the cavity {ipA^R) — 0). 

As we have shown in Ref.|6|, the time reversion for a wave packet arriving at the cavity 
from the outer region is performed through the injection 

XB{e) = ^ B , B {e)]- l ^{e). (5) 

Notice that in the above expression, 

/oo 
dtexp(iet/h)^(t), (6) 
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accounts for the time reversed complete evolution of the wave packet, i.e. both the incoming 
and outgoing components are necessary in order to compute tp B v (e). However, in the present 
work we want to reverse the signal that is produced inside the cavity. Then, we should 
deal with the building up of the complete evolution from the knowledge of only a partial 
component ^s(t) at times < t < t B . By supposing that the injection function is known, 
the reversed propagation at the boundary should be 

^B v (t) = r B (tR -t) = ih [ G BtB (t - t')x B (t'W, 0<t<t R . (7) 

Jo 

Since the PIF prescription achieves for exact reversion inside the cavity, the evolution at 
subsequent times can be assumed as the wave packet that starts at the focalization time 
with ipAitn) = ^4(0). Thus, for subsequent times we have 

4 cv (t) = ihG B , A (t - t R )fa(0), t > t R . (8) 

A comparison between Eqs.([2]) and (IE|) shows that for the most simple case in which all 
coefficients of V'a(O) are real, both evolutions become identical. The same evolution but sign 
changed should be obtained when V'a(O) is purely imaginary. 



In a more general situation where the coefficients present relative phases we start with 
the Fourier transform 
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(t-t R )r A (0)e^ h dt 
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(9a) 
(9b) 

(9c) 



where the second term within the brackets can be interpreted as the unknown evolution 
for subsequent times. At this point, we account for the hopping matrices Vb,a' an d ^a',b 
(A' C A) connecting the cavity with the outer region. The e argument in Gb,a{£) is 
neglected for a compact notation henceforth. Hence, it is possible to rewrite Gb,a in terms 
of the Dyson equation 
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and its transpose 
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where G^ -* follows Eq.(j4j) for the unperturbed Hamiltonian with Yb,a> — ^A',b — 0. For 
the considered Green's functions we have G B A = and since the subsystem A is closed and 
non absorbing, Vb,a'G a , a is a real matrix. Therefore, we obtain for the advanced Green's 



function 



and 



j b,a 



"a,b — ^a,a 



,Ya',bG* b b — Ga,b [Gb,b] [Gb,b] 



ihG B ,Aip A (0) = |-i^(0)G^ jB ] =[-ih$i(0)G At B[GB,B] '{Gb^Y 
-Gb 1 b[[Gb i b]~ 1 ]VbOO- 
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Finally, we use this last expression and rewrite Eq.( l9c|) as 
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I-G b ,b[[Gb,b} ] V>s(e) 
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Now we have the complete ip B v (e), the injection function according to the PIF formula is 
given by 

Xb{s) = ^ [GB.BpVfl'te) = ^exp(i£t R /h) Im ([G^sp 1 ) j* B (e). (15) 



This is the PIF prescription for the case where the initial state is an excitation inside the 
cavity. As we shall see, the imaginary component of the inverse in Gb,b is closely related 
with the group velocity of the scattered waves in the outer region. Therefore, contrary with 
the external source condition presented in Ref. 6j, the correction imposed by the PIF in 
internal source condition do not depends on the structure of the cavity. 

A. Application examples 

As a first attempt to evaluate Eq. (|15p we propose a one-dimensional discrete system 
whose Hamiltonian writes 

H = J2 E * \z) (x\ + ^(K-.z+i \x + 1) (x\ + h.c), (16) 

a;>0 x>0 

where E x denotes the site energy at x > and V X)X+ i is the hopping term between sites x 
and x + 1. In this system, the cavity can be delimited by only one transducer at site x$- 
The Green's function is obtained through the continued fraction technique as 

£ — E xs — Li n (e) — L ou t(e) 
where E in (e) and E out (e) are the self-energy corrections to the site energy E xs due to the 



presence of the remaining sites inside and outside the cavity respectively 111 ]. 

Decimation on the cavity gives the T,- m (e) contribution as a continued fraction whose 
parameters are the hoppings and the site energies and, in absence of magnetic fields and 
absorption processes, they are all real. Therefore, 

^in(£) = Vx s ,xs-1 I ^x S -i,x s , 

e ~ ^x S -i ~ Vxs-\,x S -2 1 V Xs _ 2)Xs _ 1 

(18) 



£ — -c/n 



is also a real function. 

In the thermodynamic limit where the number of sites outside increases indefinitely, the 
homogeneous outer region gives the self energy as a complex number: 

E out (e) = A <rat (e)-iT out (e), (19) 



with 

' ^ + \J{^T-f ~ V ' 2 e-E Q <-2 \V\ 



A out (e) = A(e-E Q ) = < 
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\e-E \ <2\V\ , (20) 



e-Ep _ a / f £ -Ep \ 2 



(?=£*) -V* e-E >2\V\ 



and 



^ = r (£ - Bo ) H ^-^) l«-*>l^, (21) 

|e-£ |>2|V| 

whose parameters V and £?o are the hopping and the site energy in the outer region respec- 
tively. The relationship between r out (e) and the group velocity can be found through the 
dispersion relation for the asymptotic waves: 

e k = E - 2Vcos(ka), (22a) 



v ° = n -dk = x asm(A;a) = TV - {—2—) (22b) 

= ^r out (e fc ), (22c) 

n 

with a the distance between first neighbor sites. Remember that choosing Eq = 2V so that 
the differential form of the Schrodinger equation for a particle with mass m is obtained as a 
limit h 2 /2Va 2 — ► m when V — > oo and a 2 — » 0. Since the only imaginary component from 
Sout(^) is required, Eq. (IT51) writes 

X pw(x S , e) = X b{x s , e) = ^e i£t ^ h T out (e)r(x s , e). (23) 

The complex exponential serves to define the origin of time, and hence we can neglect it. 
Using the definition of Tie) given above we get 



X pi F (x 5 ,5) = t/1 - ( 1 ^) 2 f nx 8 ,e) (24a) 



e - 2F X ' 
2V 



Xtrm(zs,£)> (24b) 



perfect inverse filter 

i.e., the time reversal mirror prescription of injecting an amplitude proportional to the 
detected wave function must be corrected by a filter that enables the exact reversal. When 
the initial state is a local excitation composed by a few sites inside the cavity, the Fourier 
transform of the detected signal i[}(xs,e) will cover the whole energy band between — 2V 



and 2V. The factor in the PIF procedure shows that the correction becomes effective near 
the band edges where r out ~ 0. However, since in simple homogeneous cases the filter is 
just proportional to the group velocity, it is clear that in the acoustical case Xpif(%s, £ ) — 

Xtrm(xs,z)- 

A more general case where Eq. fll5j) can be tested is the multichannel system schematically 
shown in FigHJ Here, a two dimensional cavity is surrounded by Nl transducers at the left 
and Nr at right. The total Hamiltonian is written as 

( M L V L (e) N 

M c V fl (e) , W) 

V£(e) M R 

where He is the Hamiltonian inside the cavity and Ml (Mr) is the Hamiltonian on the left 
(right) side. Vc(£) and Yr(e) are real matrices containing the hopping terms between the 
boundaries and the cavity. The unitary transform 



V£(e) 



V 



/ 



/ 



u 



V L 
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V 



(26) 
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l c o 
V R 

with lc the identity matrix in the cavity, separates the channels into Nl + Nr indepen- 
dent one-dimensional homogeneous chains. The energies E% follow the dispersion relation 
Eq.f l22ap for the discrete case: 

/C7T 



E* = Ec 



2Vcos 
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U nk 



U kn 



l,...,N a . 



nkir 



sin 



(27) 



(28) 



N a + 1 \N a + l, 
By performing a decimation process in each one of these chains, the effective site energy at 

the boundaries is 

[UlH Q U a ] kkl = [Jgff - E(e - B%)\ 5 Kkl) (29) 

where a = L,R and k,k' G B denotes the eigenstate index in the transducers subsystem. 
We use again the decimation procedure on the cavity sites in order to obtain the total 
Hamiltonian in the reduced basis of the boundaries Hss(e). Since the Green's function is 
defined as: 

[(el - My 1 ] B B = U (eI B , B - H b ,b(£)) -1 U f , (30) 



*B,B 
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we have 



[G^b] -1 = U [eI B , B - M BiB (e)} U f , (31a) 

Im [G BiB ] _1 = -Ulm [H B)B (e)] U f . (31b) 

The only imaginary components in M. BtB (e) come from the diagonal part and, as indicated 
in Eq.§SD, 

Im [G B , B ]^ l)ea = £ < fc r ( £ - E kKi, (32) 

k 

where (i,j) G a means that both n and I are taken in the same side of the system: mixing 
terms between L and R sides of the above matrix should be zero. Thus, the injection 
function is obtained replacing Eg. (I3"2"|) in Eg. (TT5"l) as, 

2 Na 

[XB, a (£)l = jf etR ' h E <^ £ - E %Xi &,«( e ) . • ( 33 ) 

k,l=l 

As we have seen previously, the PIF correction for the internal source case is given by 
the group velocity at the boundaries. Notice that the r(e) factors are centered on the eigen- 
energies E% whose width is given by the hoppings in the outer region. As we show in the 
numerical simulation, depending on the particular choise of these parameters the correction 
would be non trivial. 

B. Numerical simulation 

In the multichannel system depicted in Fig{T]we want to reverse the dynamic that follows 
an initial condition of local excitation in the site labeled as Xq. We perform a numerical 



simulation for the time domain using the Trotter-Suzuki algorithm 12j. In such case we 
have Nl = 3 transducers for the left and N B = 2 for the right and the propagation is favored 
through the vertical hoppings (according to Fig{T]) of the channels whose values are twice 
the horizontal ones. 

The time-dependent response functions G(xs', t;xs, 0) between the transducers at the 
boundaries are then calculated and a numerical Fourier Transform is performed in order to 
obtain the 5x5 matrix G BjB (e). A comparison between analytical and numerical evaluation 
of Eg. (1321) is shown in Fig f5] for some representative elements. They show how the group 
velocities in each channel are mixed. Figures [2^, [2b and [2b show the diagonal elements and 



the Fourier transform of the detected signals for transducers 1,2 and 4 respectively. For 
transducers 2 and 4, the correction becomes effective in the central region of the spectrum 
where the detected signals are appreciable. As we observed from Eg. (152"]) . Figf2fe shows that 
no mixing exist between transducers that belongs to different sides of the cavity. We can 
observe a little portion of points from the numerical results that does not fit with the ana- 
lytical curves. These errors are related with the approximation done by the Trotter-Suzuki 
and because a finite registration time was taken. Once we have the corrected injection, we 
reverse the dynamic and compare it with the TRM case. Although TRM achieves a reason- 
able quality in the spatial focalization, differences appear both for the time evolution at the 
focal site (see FigJH]) and the spatial distribution of the recovered signal at the focalization 
time tft shown in FigJU Notice that in the case where the eigen-energies are close together 
and the band width becomes wider, the PIF correction is only a constant pre-factor and the 
focalization quality by TRM increases. 

III. COMPARISON WITH OTHER RESULTS 

In a more general point of view, we may develop an alternative derivation of the PIF for- 



mula based on the fundamental optical theorem [13( . In words, it states that the continuous 
density of states at a given site x is built upon the escape rates through the boundaries of 
the system. If we chose a semi-infinite chain, there is no escape to the left of site xq while 
an excitation can escape to the right of site x$ with a rate Im.Il xs /h. In equations: 

Im Gx ,x = G X0:XS Im E Xs G xs>XQ . (34) 

Multiplying both sides by a pulse s(xo,s) and the ih factor one finds 

ifr (G XOtXo - G* XQ>X0 ) s(x ,e) = 2iG X0tXS T xs (-ih)G* XS;XQ s(x ,e). (35) 

If we identify the outgoing wave originated on the excitation at the focal point as 

ip(x , e) = ihG XOtXo s(x , e), (36) 

we re-obtain the PIF formula by identifying the second term in the left hand side as the 
incoming wave function that produced the pulse in the focal point. The complete evolution, 
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which is equal to the time reversed one, results 



ip(x , e) + ^*(x , e) = Aev( x o, e ) (37a) 

s 1 



l"-tja;o,a;c 



-ra, 5 V*(zs,e) 



(37b) 



Xpif 



The term in square brackets is the excitation that must be injected at xs to produce -^rev- 
Notice also that this prescription does not depend on x . 

It is interesting to show some similarities between PIF and the focusing procedures (STIF 
and miSTIF) described in Ref.[lOj. We will consider a simple example with a single trans- 
ducer at xs and a focusing point located in x . The goal of STIF is to inject some signal at 
Xs that, at a later time, is able to form a spatio-temporal delta function in xq. The strategy 
of STIF is to compute the inverse of the propagator G(xs,xo;e) = G XStX0 . The injection 
proposed by the STIF could be interpreted as that which achieves the time reversal of an 
hypothetical pulse emitted by a source at x , i.e. s(x, t) = S x>Xo S(t), resulting 

Xstw(x s ,s) = i^G~ 0iXsW s(x ,e) = — [G* XStXO G XOtXs ] G* XSjXoW s(x ,e). (38) 

Since the time reversal mirror injection is given by 

Xtrm(x s ,£) = ^*(x s ,e), (39) 

and 

^*(x 5 , e) = -ihG* XStXQ s{x , e), (40) 

we can rewrite the STIF as a filter 

Xstif(x s , e) = -^ [G* XS)XQ G XOiXS ] ~ l Xtrm^s, e). (41) 

Notice that this filter requires the determination of the non-local Green's function connecting 
the transducer at xs and the focal point xq. Therefore, the miSTIF was developed as a 
procedure that replaces it in terms of the local Green's function at the transducer site. Such 
replacement is enabled by the optical theorem 

Im G XStXS = G XSjX0 Im ^ Xo G* XOjXs + G XSjXS Im S a , g G* s X5 , (42) 

applied to an infinite chain where xs > Xq. Here, S~ o (£1 „) is the self-energy correction that 
accounts for the free propagation at the left (right) of the focal point (transducer) whose 
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imaginary part is V XQ (r+ g ). Since 

lmG XSjXs = Gx s ,x s W^atf + ^rrJG^s,^ (43) 

we may use that, in a homogeneous system, the escape at both sides has equal group 
velocities 

Im[£+ s + £; s ] = -2IV (44) 

In consequence, under these restricted condition, the optical theorem leads to 

from which the miSTIF prescription can be obtained as a broad-band approximation to 
Eq.flHD in which T xo (e) = T Xs (e) ~ T = V, 

XmiSTl?{xS,£) = 2V [ - ^ G x a *s\~ X Xtrm^s, e), (46) 

i.e. the minimally invasive spatio-temporal inverse filter. 

A comparison between miSTIF and PIF shows that both coincide since the self-energy 
correction to the left and right sides of the chain are the same and 

Gxs ' xs = e-2A(e) + i2T(e) (47a) 

(47b) 



2r(e)' 

and hence, in that particular case 

[ImG^r^-Im^J. (48) 

and we see that both injection prescriptions, the PIF and miSTIF, are proportional. 

In turn, since PIF is an exact procedure, it allows to test the efficiency of miSTIF in 
focusing the excitation to achieve a spatio-temporal delta function. This actually results 
as the back-propagation of the "original" pulse excitation. As such, in our one dimensional 
homogenous case, the amplitude at the focal point is not a delta function in the time domain 
but Jo(2Vt/h), a Bessel function of first kind. In the acoustic case, where the wave function 
normalization is not required, this can becomes an approximation to the delta function as 
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IV. CONCLUSIONS 

We presented the generalization of the PIF focusing method in the internal source con- 
dition. Although we considered the Schrodinger equation, since this formulation is based 
in general properties of the propagators, as before 8|, the results remain valid to describe 
sound propagation. Contrary to what happens when the source is external, the prescription 
one has to use in this case does not involve internal details of the cavity but only simpler 
information about the propagation in the outer region. This simplified filter applies to two 
physically relevant situations: 1) When the excitation is actually originated in the interior 
of the cavity by an external source. In this case, one obtains a perfect reversal of the wave 
function for all the times after the source have been turned-off. 2) Whenever one uses an 
external source, but one can separate the incoming from the outgoing waves. In that case, 
one can use the recording of the outgoing wave to perfectly reverse the whole excursion 
of the excitation through the cavity. In principle, one can achieve this condition when the 
boundaries are placed far enough from the reverberant region. An interesting consequence of 
our result is the affirmation that for close cavities in the acoustic case, TRM should produce 
perfect time reversal, up to a constant factor, provided that the external propagation is free 
of collisions. 

We confirm the validity of this expression both analytically and numerically in a multi- 
channel system and compared with its TRM counterpart. Results on the reversion quality 
show that PIF correction constitute a notable improvement over TRM prescription in cases 
where the energy dependence has a non-linear dispersion relation, as typical excitations de- 
scribed by the Schrodinger equation, or when the escape velocity presents structures due 
to the presence of many simultaneously propagating channels or when collisions outside the 
cavity have a relevant contribution. 
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FIG. 2: (Color Online) Representative elements of the PIF correcion: numerical in black points 
and analytical in blue line. Real and imaginary components of ?p(xs,£) are shown for comparison 



reasons. 
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FIG. 3: (Color online) Recovering of the initial condition in the site xq. Black line: ideal reversed 
evolution. Red line: normalized TRM. Blue line: PIF. 
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FIG. 4: (Color online) Normalized spatial distribution for the recovered initial condition in the 
cavity. The location of the transducers is indicated by rectangles. 
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